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Abstract 

o 

f"^ ■ We use the 2-dimensional a model as a toy model to study the behavior 

of anomalous amplitudes in the limit where the constituent quark mass is 
small. Symmetry arguments tell that the 7r° — > 7 amplitude should vanish 
if m — * 0, but we show that this conclusion is spoiled by infrared singu- 
|_ larities. When a proper regularization (resummation of a thermal mass, 

for instance) is taken into account, this amplitude vanishes as expected. 
We also study the amplitude n°a — > 7 and show that it does not vanish 
in the same limit. 
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1 Introduction 



In the past years, many works have been devoted to the study of anomalous 
processes at finite temperature and density, and in particular near the chiral 
symmetry restoration. In particular, Pisarski jj], ^ concluded that the neutral 
pion decay amplitude ir° — + 2j vanish when chiral symmetry is restored^]. This 
conclusion is based on a direct calculation of the corresponding diagrams at 
finite temperature (and zero density) in the imaginary time formalism, and on 
symmetry considerations that forbid certain couplings in the symmetric phase. 
An additional conclusion was that this decay might be replaced by Tr°a — > 27 
in the chirally symmetric phase (which is allowed by the same symmetry argu- 
ment), as indicated by a calculation of the box diagram at finite temperature. 

This result has been confirmed by Baier, Dirks and Kober (jllj and by Salcedo 
| fl2"| using functional approaches in which one integrates out the fermions at 
the level of the the generating functional. Salcedo in 112] gives also general 
arguments according to which 7r° — > 27 (or 7r° — > 7 in two dimensions) should 
vanish in the chiral phase, and be replaced by amplitudes involving the a meson. 
From a technical perspective, the common point of all these studies is the use at 
some stage of the imaginary time formalism, in the limit of vanishing external 
momenta. 

In another study, one of us |k| studied how the neutral pion decay amplitude 
depends on the kinematical configuration of the external legs in order to explain 
the discrepancies found between |^] and calculations performed in the real 
time formalism by @ [u|. To that purpose, the 7r°77 amplitude has been 
calculated at finite T in the real time formalism, in the limit of small external 
momenta. It appeared that this limit cannot be uniquely defined (it depends 
on the path followed to reach the zero momenta point) and that the results 
of pi ^, pi [l^] concerning this amplitude do not correspond to its on-shell 
value, but to a different way of reaching the limitQ. Additionally, the conclusion 
according to which the pion decay amplitude vanish above the critical point 
appeared to be questionable since the physical (on-shell) amplitude has a non 
vanishing limit at the critical point. 

To accommodate this result with the general arguments provided in |], ^ |i2f , 
one can notice that both Pisarski's symmetry argument jl], ^ and Salcedo's 
argument |T^| amount to the fact that one power of the quark mass (the mass 
the quarks acquire through the spontaneous breakdown of chiral symmetry, via 
their coupling to the average value of the a field) appears in the numerator when 
evaluating the Dirac trace associated to anomalous amplitudes. Therefore, since 

1 This does not contradict the wdl established fact that the coefficient of the axial anomaly 
is temperature independent (see M, UHH| for instance). Indeed, it has been shown in 

§^j, |l(j] that the relationship between the anomaly and the amplitude can be modified by 
e existence of an additional 4-vector (the 4-velocity of the plasma in the observer's 
frame) that can enter the general form of thermal amplitudes. 

2 The reason for this is easy to understand: since the energy variables are discrete in the 
imaginary time formalism, the only way one can consider the "zero momenta limit" in this 
formalism is to set first the discrete bosonic energies to zero, and then take the limit of zero 
three momenta. This way, the external momenta are forced to be space-like. 
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the average value (a) goes to zero in the chiral phase, the numerator vanishes 
above the critical point. Implicit in the argument is the fact that the correct 
dimension is provided by inverse powers of the temperature (as opposed to 
powers of the quark mass), as is the case in the imaginary time formalism at 
the static point. In other words, this argument is valid only if the denominator 
does not vanish when the mass goes to zero. This is precisely what fails when the 
amplitude is calculated on shell. One gets as expected one power of m = g{a) 
in the numerator, but the denominator turns out to be 1/mT. 

The remaining power of m in the denominator indicates that the infrared or 
collinear behavior of the triangle diagram worsens when m — ► 0. In fact, as noted 
in jl3| and fl^j , the constituent quark mass m ceases to be the relevant infrared 
regulator when m is smaller than gT, and should be replaced by a thermal mass 
of order gT that do not vanish in the chiral limit. Since an additional property of 
fcrmionic thermal masses is that they respect chiral symmetry, this thermal mass 
cannot appear in the Dirac's trace. As a consequence the result m/mT obtained 
for the on-shell amplitude in the bare theory becomes m/mthT after one has 
resummed the quark thermal mass mth ~ gT (if m <C mth)- The consequence 
of this regularization is that the resummed on-shell decay amplitude vanish in 
the chiral limit. In other words, Pisarski's symmetry arguments holds for the 
physical amplitude only after a proper regularization (in order to get rid of all 
potential infrared or collinear singularities) has been issued. 

Essential in this discussion is the influence of the kincmatical conditions 
for the external legs on the infrared behavior of a thermal amplitude, since it 
can dramatically alter one's conclusions. The second important point is that a 
calculation in the imaginary time formalism at the static point does not give a 
physical amplitude. Therefore, it would be interesting to test the second half 
of Pisarski's conclusions, related to the ir°a — > 2j amplitude, by calculating 
this amplitude in the real time formalism and studying how it depends on the 
kinematics (up to now, this amplitude has only be considered at the static point). 
Since this amplitude is given by a four point function, it is a very complicated 
task to extract this behavior in its full generality. There is however a toy model 
in which this kind of study can be done quite simply: the 2-dimensional linear a 
model. Indeed, in this model the neutral pion decays into a single photon, and 
the analogous of the 4-point amplitude suggested by Pisarski would be ir°<j — ► 7, 
which, being a 3-point function, is rather easy to calculate. 

The present paper is devoted to an analysis of the anomalous amplitudes in 
the 2-dimcnsional a model at finite temperature and chemical potential. We 
consider both the pion decay ir° — > 7 2-point function, and the n°a — > 7 3- 
point function. Emphasis is put on studying how these functions depend on 
the kinematics in the limit of small external momenta, near the chiral limit (m 
small compared to /1 and T). 

The structure of the paper is as follows. Section || defines the model, as 
well as some notations and shorthands that will be used extensively later. In 
section |[ we calculate the amplitude for the n° — ► 7 decay, and reduce it 
to a very compact form. In section ||, we study the amplitude of tt°<7 — ► 7. 
Although a priori much more involved, this amplitude can also be reduced to a 
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simple expression. All our results are expressed in terms of some function I(K) 
defined by an integral. Basic properties and limits of this function are derived in 
appendix |A|. Finally, appendix |^ derives some relations between a few integrals 
that appear in intermediate stages of section 

2 Conventions and notations 

We consider the 2-dimensional linear a model [[l8[ with two quark flavors, 
in which the mesons are coupled to quark fields as indicated by the following 
Lagrangian: 

£ = i^]py -2g^(at +iir-t"f 5 )^ , (1) 

where t = 1/2 and Tr(t a t b ) = S ab /2. We recall that in two dimensions the 
Dirac algebra is defined by the following set of relations: 

7 5 - \e» v Wu , (2) 

where e^ u is the 2-dimensional Levi-Civita tensor, normalized by e 01 = +1. For 
later use, let us quote first a generic trace formula: 

Tr(#$7V) = A • BTr«? 7 V) - A ■ CTr($ 7 V) + B ■ CTr(4 7 V) . (3) 

In order to keep the following expressions compact, it is helpful to define the 
"dual" of a given vector by: 

A" = <? V A V , (4) 
as well as the "wedge product" of two vectors: 

Af\B = eTApBv . (5) 

According to these definitions, we have the obvious relations 

A = A 

AAB = A-B 

(A AS) 2 = (A-Bf - A 2 B 2 . (6) 
When A + B + C = 0, we have also: 

AAB = B AC = C AA . (7) 

Finally, we have 

Tr(4 7 V) = -2A^ . (8) 
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Figure 1: One- loop ir°j amplitude. 

3 7T° — ► 7 amplitude 

3.1 Retarded amplitude 

We consider first the 1-loop contribution to the ir° — + 7 decay amplitude de- 
picted on figure [l| The Feynman rules for the retarded-advanced formalism^ 
[ [20I pl| p2| give for the retarded amplitude the following expression 

n»(K) = -eg J ^Tr((^ + m)7^(^ + ^ + m) 7 5 ) 
x{n F (/ ,M)DiscA R (i) A R (i + if) 
+n F (l + k Q ,n)~DiscA R (L + K) A A (L)} , (9) 

where A H A (L) = i/(-^ 2 — to 2 ±iZoe) are the retarded and advanced propagators, 
n F (lo, n) = l/(exp((?o — /i) /T) + 1) is the Fcrmi-Dirac distribution function, and 
where the notation "Disc" denotes the discontinuity across the real energy axis: 

DiscA,, (L) = A R (L) - A A (L) = 2^e(Z )<5(L 2 - to 2 ) . (10) 

The expression of the trace is very simple 

Tr ((£ + m)^{t + + m) 7 5 ) = -2mK» , (11) 

and in particular it makes obvious the fact that II M is transverse with respect 
to the photon momentum. Moreover, being independent of the loop momentum 
L, it can be immediately factorized out of the integral. 

3.2 Zero momentum limit 

At this point, it is convenient to perform the change of variable L + K — > — L 
on the second term of Eq. (^) in order to make the expression more symmetric. 

3 One should pay special attention to the chemical potential. Indeed, one should use a 
chemical potential — /i in statistical weights where the Feynman rules give an argument — feo- 
In other words, /i appears in the formalism to account for the fact that the fermions carry 
some conserved charge, and the sign of this charge for a given propagator depends on how 
one orientates the propagator. 
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Then, the Dirac distributions hidden in the discontinuities make the integration 
over Iq trivial, which gives^J 

Km _ _ w 7^ "' ( -"-t"' ( "-"" £ ■ <"> 



where we denote = V^ 2 + m 2 and = (ryw;, /). It is now trivial to perform 
an expansion in powers of the external momentum K. The first term in this 
expansion, of degree in K, vanish because the corresponding integrand is an 
odd function of I. The first non vanishing term in this expansion comes at the 
next order, and is of degree 1 in K: 

T1%(K) = imegK^I{K) , (13) 

where I(K) is an homogeneous function of degree in K , containing the non 
trivial part of the momentum dependence, and defined by: 

I{K) = J 2^ (ZTW ■ (W) 

— oo 

This integral is studied in some important limits in the appendix |X| 
3.3 Discussion 

We observe for this amplitude the same features as in four dimensions. The 
most striking effect is related to what happens near the chiral limit m — ► 0, 
and is visible in formulas Eqs. j3l| ) and (|38|). In the limit m/T, mj \i — ► 0, the 
function I(K) behaves as follows: 



k 




8mT ViPcosh 2 (/x/2T) 



if m<T 

If fc = , I(K) = { * . (15) 

if T</i 



4717* 2 



The configuration obtained with fco = and /i <C T corresponds to Salcedo's 
result, previously obtained in the imaginary time formalism at the static point. 
This point is particular because the first term in the expansion of I(K) at small 
m vanishes if fco = 0. For any other point, the expansion starts with a term 
behaving like I(K) ~ 1/mT. However, as one turns the chemical potential on, 

4 We have dropped the R/A prescription for the denominator, since it can easily be recov- 
ered at the very end of the calculation by substituting fco — * kg + i0 + . 
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we see that this leading term is exponentially suppressed when fj, 3> T. As a 
consequence, in a dense and cold system, the function I(K) starts by a term in 
l//x 2 , whether ko = or not. 

The reason why the on-shell value is so singular when m — > is related to 
collinear singularities: in 1 + 1 dimensions, all the spatial vectors are aligned, so 
that we are always at the most singular point 5 . Here also, in order to be able to 
apply Pisarski's argument, one should first regularize the theory by resumming 
a thermal mass[jl3[ [l7]]. Then, all the powers of m in the denominators would 
be replaced by powers of the thermal mass, leaving an uncompensated power of 
m — > in the numerator. From Eq. (|l3|), we can write an effective Lagrangian 
coupling the neutral pion to the photon: 



C o 



egme^ J d 2 xA"(x)I(id x ) d v x -K°\x) 
eg\ u I ' d 2 xA»(x)I(id x ) (a)d^°(x) 



(16) 



This result completes the effective coupling found by Salcedo in [12[] , by incorpo- 
rating all the nonlocal terms. The reason why the non-locality of this coupling 
has been missed in Jl2[ can be traced back in a misuse of the imaginary time 
techniques to get the zero momenta limit. 

4 7T°a —> 7 amplitude 
4.1 Retarded amplitude 

We now consider the one-loop contribution to the ir a — > 7 amplitude, repre- 
sented on figure 0. In the retarded-advanced formalism, the component 




0\ o 



+ Lv 








Figure 2: One-loop ir°a^ amplitude, 
of the vertex receives the following contribution from the first diagram: 

T^(K, P, S) = ie .g 2 J -^^{{t + m)7 5 (^ - P + m)<f& +$ + m)) 

5 This is to be contrasted with what happens in four dimensions: there the collinear singu- 
larities are softened by subsequent angular integrations, so that the on-shell amplitude is not 
exceptionally singular JlSj, 16|. 
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x [n F (Jo, /i)DiscA a (L) A R (L + S)A A (L - P) 
+n F (l Q - po, M)DiscA R (L - P) A R (L)A R (L + S) 
+n F (l + so,n)T>iscA R (L + S) A A (L)A A (L - P)} . (17) 

Again, the expression becomes simpler if we perform the changes of variables 
L — P — > L on the second term, and L + S —* L on the third one. This enables 
one to have common statistical weight and discontinuity for the three terms, 
the tradeoff being that the trace becomes different for the three terms. We can 
apply the same manipulations to the contribution of the second diagram. With 
the additional change L — > — L on the second diagram, we can merge the two 
contributions and find 

+ OC 



{ 



?7=±1 

Tr£ Tr£ 



(2L V ■ K + K 2 )(-2L n -S + S 2 ) {2L n ■ P + P 2 )(-2L n ■ K + K 2 ) 

+ (2L V -S + S 2 )(-2L 71 -P + P 2 )i ' (18) 
where the traces are given by 

Tr£ = -2m 2 TrGP 7 V) - 2L n ■ STr^V) - Tr(A^ 7 V) 

Tr£ = -2m 2 Tr(/> 7 V) + 2L, • PTr(A, 7 V) - Tr^/^V) 

Tr£ = -2m 2 Tr(/> 7 V) - Tr^/W) . (19) 

4.2 Zero momenta limit 

We can now proceed with the expansion in powers of the external momenta. The 
first term in this expansion is of degree —1 in the external momenta. An explicit 
calculation of this term shows that it vanishes thanks to energy-momentum 
conservation: P + K + S = 0. The next term, of degree in the external 
momenta, vanishes also because the corresponding integrand is an odd function 
of I. Therefore, the first non vanishing term is of degree 1 in the external 
momenta. An explicit extraction of this term gives 

T %R (K,P,S) = 4 / ^ M^i^ E { _ m2Tr( p 7V) 

— oo 77— ±1 

S 2 S 2 K 2 K 2 P 2 



(L v ■ P) 2 (L v ■ S) 2 (L v ■ S) 2 (L v ■ K) 2 (L v ■ K) 2 (L v ■ P) 2 
1 / K 4 P 4 S 4 



+ 



(V ^)(V P )(V 5 ) V V ^ L v ■ P L v 



-) 
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K 2 



S' 4 



+ 



L v -K Vl v -K\L v -P L v -SJ ' {L v -S) 2 {L v -P} 2 
Tr(A^ 7 V) r S 2 K 2 



(L v ■ K)(L ri ■ S) 
Trifcp^j 5 ^) r K 2 P 2 



L^i ' K L/fj ■ P 



(L v ■ P)(L n ■ K) 



(L v ■ S)(L V ■ P) 



} 



(20) 



At this point, we can make use of Eqs. (|) and (|). It is now obvious that the 
result can be expressed in terms of the following integrals^ 



+ oo 



dl n F (-uj h ^) -n F (uj h fj 1 ) 1 1 
2?r 2ui L+- AL+-B 



m = 

AAB — 



dl n F (-uii,fi) -n F (uii,fi) 1 _ I{A) 
2ir 2uji (L+ ■ A) 2 ~ A 2 

O 

dl n F (- uji, n) - n F (o>;, /x) i+ 1 
2tt 



2wj 



(L+ ■ A) 2 L + ■ B 



-|-oc 



= m 2 



dl n F {-wun)-n F {ui,n) 1 



1 1 



2tt 



2u>i 



{L+ -A) 2 L + -BL + -C 



-\-oc 



J, 



= m 2 



dl n F (-wi, n) - n F (u>i,fj,) 1 



2tt 2w; 

— OO 

where A, _B, C can be any of P, 5, K. 



(L+ ■ A) 2 (L+ ■ B) 2 



(21) 



4.3 Trans versality 

The major difference compared to the case of the 7r°7 amplitude is that the 
above three-point amplitude is not manifestly transverse with respect to the 
photon momentum. The fact that the Dirac's traces depend upon the loop 
momentum L indicates that it could be necessary to calculate all the integrals 
before one can see the transversality. The situation is not so intricate though, 
since it happens that we only need to establish some relations between the 
various integrals defined in Eqs. (f2l|). In appendix we show how the last 
three integrals of Eq. can be expressed as functions of the first two. Using 
the relation given in Eqs. (f4(i|), ( [43|) and as well as Eq. (ff5|), it is a simple 

6 Note that we need only J AABC when A + B + C = 0. Having calculated J AABC under 
this restrictive assumption, one cannot obtain J AABB from it by enforcing C = B. 
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matter of algebra to check the transversality of the 7r°<77 amplitude with respect 
to K, without the need of calculating the various J AA and J AB . 

Once we know that the result can indeed be written as r M = YK^ : one can 
contract the amplitude with S for instance in order to extract the coefficient T. 
A straightforward calculation gives 



with 



T» RR (K, P, S) = teg 2 K»F(P, S) , (22) 

ev d cn - P 2 S 2 I(S) + (P ■ S)[P 2 I(P) + (P ■ K)I(K)} 
F(P, S) = , (23) 

where implicitly K = —P — S. At first sight, this expression could explode 
whenever two momenta become parallel. However, one can check that this is 
not the case, because the numerator behaves like (PAS) 2 when PAS becomes 
small. 

4.4 Discussion 

We see that this 3-point amplitude involving the a field depends on the same 
function /(.) defined above, and has one power of the mass m less when com- 
pared to the 7r°7 amplitude, in agreement with the general arguments of [jj], ^|. 

Again, it is found that this limit depends on the kinematics, i.e. on the way 
one is approaching the zero momenta limit. In particular, the way this amplitude 
depends on m at small m depends strongly on the kinematics. This is to be 
contrasted with the result of p^ |, which only picked one particular limit. Except 
at the static point (fco — po = sq = 0), this amplitude becomes singular at the 
critical point where m — > 0, indicating the necessity of regularizing the fermion 
propagator by a thermal mass. After this resummation has been performed, 
this amplitude is regular but does not vanish in the chirally symmetric phase. 
Therefore, the conjecture of H holds, but only after infrared regularization. 

One can also write an effective coupling associated with this amplitude: 



£„o ay =eg e M „ / d xA fl (x)F(id Xl ,id X2 ) 

x[a(x 2 )d:y(x 1 ) + 7T (x 1 )d: 2 a(x 2 )} . (24) 



£1— £2— £ 



Of course, if one uses Eq. (|25|), one finds a local limit for this effective coupling 
in the limit of zero temperature and density. 

5 Conclusions 

In this paper, we have studied the 7r°7 and it°a~/ amplitudes in the 2-dimcnsional 
a model at finite temperature and density. For both of these amplitudes, the 
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zero momenta limit is not unique and strongly depends on the kinematical con- 
figuration. In particular, the imaginary time formalism should be used with 
great care when looking at this limit. Indeed, if one first sets the external dis- 
crete energies to zero, then all the information regarding the non-locality of the 
amplitude is lost, and in particular the physical limit cannot be recovered. A 
proper way to use the imaginary time formalism would be to perform the sum 
over the loop discrete energies while keeping nonzero external discrete ener- 
gies. After that, one should perform the analytical continuation to real external 
energies, and only then consider the zero momenta limit. 

The other conclusion of this work is that collinear or infrared singularities 
spoil the general symmetry arguments given in ||, ||, [l2| to justify the nullity of 
the pion decay into photons in the chirally symmetric phase: the on-shell decay 
amplitude in the bare theory does not vanish. For these arguments to be valid, 
one should perform the resummation of a thermal mass that will regularize the 
fermion propagators. 

If such a regularization is used, then the conclusion is that 7r° — > 7 vanishes 
when m — > 0, while ir°a — > 7 does not, in agreement with the conjecture of 
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A Properties of the function I(K) 
A.l Vacuum limit 

The purpose of this appendix is to study the integral I(K) since all the quantities 
calculated in this paper can be expressed in terms of this function. A first check 
is to look at the zero temperature and chemical potential limit of this function, 
which gives immediately 

lim I(K) = —L, . (25) 
7>--o+ 27rm 2 

As one could expect, this T = fi = limit does not exhibit any non-locality in 
its momentum dependence, since this is a purely thermal feature. 

A. 2 Transformation into a sum 

A convenient way to look at the high temperature or density limit is to turn the 
integral defining I{K) in Eq. Jl^) into a sum, by making use of the following 
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identity^: 



1 1 +^ 

+ 1 2 ^ 



x 2 + (2n + 1) 2 tt 2 

n— 



(27) 



This identity enables one to rewrite n F (— u>i,fi) — n F (u>i,fi) as a series, from 
which it is straightforward to first obtain 



I(K) 



1 K' 



ttT 2 k 2 



n=0 



(/(.' 



« 2 + A 2 , 



(w 2 +e)2 {v2+A 2 n) 2 + B 2 



(28) 



where we denote 



k = ko/k 



m k 



A 



m 2 jj," 



B, 



4^ 2 (2n+l) 2 ^. (29) 



At this stage, it remains to perform term by term the integration over dv, which 
is elementary and yields 



dv 



v 2 +e+ 



A 2 



<:) 



(v 2 +e-) 2 (v 2 + Al) 2 + B 2 2£_S£_ 



i (e + -e)(Ai-e) 
e_ (Ai-er+B 2 



Ai+Bi- 2e + Ai + (e - e ) VK+m + ae 



Bi 



VK + Bj + Al 
2(Ai + B 2 n ) 



(30) 



A. 3 Chiral limit far from the light cone 

Another interesting limit is the chiral limit, where m/T goes to zero, while li/T 
is kept fixed. It is very easy to extract from the above formula a systematic 
expansion in powers of m/T. It is just a matter of expanding at small £± the 
above expression, which gives for the first two orders: 



I(K) = \ 
ml 



Ik 2 ' 



m K 2 + 2fc 2 (^\, n ( m2 
T K 2 \t) XT 2 



(31) 



7 To derive this formula, one can start from Mittag-Leffler's expansion of the cot function 



cot(z) 



-t-oo 



(26) 
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where the coefficients are given by 



Fn 



a 

T 



2 +00 



— ; (2n + 1 + t^/ttT) 2 



n=0 



1 +00 



(2n+l) 2 -3^ 2 /7r 2 r 2 
((2n + l) 2 +A* 2 /7r 2 r 2 ) 3 



n+ 1) 



n=0 



Introducing the "digamma" function 



dz 



\nT(z) 



and a series representation of its first derivative 



(32) 



(33) 



1>' 



1 + z 



(2n + 1 + z) 2 

n— 7 



(34) 



it is immediate to check 



1 

8^ 



1 . I* 
h Z 

2 2ttT 



1 



32tt 3 

703) 
8tt 3 
T 2 



1 2 2ttT 



1 

2 



2ttT 



1 



Scoslr (/j/2T) 



1 ■ M 

2 Z 2ttT 



47t/j 2 



if /J<T 
if T <C /x 



(36) 



A. 4 Light cone limit 



Note however that this expansion is not valid near the light cone. Indeed, its 
derivation assumed that a small m/T would imply a small which is not true 
if K 2 is small (or equivalently n 2 ss 1). Sufficiently close to the light cone, £_ 
becomes large and a different kind of expansion must be considered. In this 
region of phase space, one can write: 



lim I(K) 



= lim — 

?_^+oo 2£ 



1 2_i^ 

_T 2 k 2 - 1 ^ ' 



1 



lim /o ^ 7 



T1 = 

2 2tt 



[£_ + (2n + lh] 2 



3 The equality on the first line is exact and comes from the formula |24| 

r(z)r(i - z) = , 

smiirz) 



(37) 



(35) 



while the limit )i ^> T in the last line is obtained from Stirling's expansion J23J ^| for V(z). 
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where we have used (k 2 — 1) 1 w £ 2 T 2 /to 2 . Making use of Stirling's formula 
[f|§ jMj, we obtain: 

lim I(K) = -l-j . (38) 

X 2 ^0 27TTO 2 

We notice that the on-shell value of I(K) is totally immune to corrections due 
to temperature or density. The origin of this property can be understood from 
Eq. (pr|): when K 2 = 0, the denominator behaves like fc 2 TO 4 /Z 4 (this is rem- 
iniscent of a collinear singularity cured by the mass m), and the integral is 
completely dominated by its ultraviolet sector. As a consequence, I(K) is sat- 
urated by the vacuum contribution for this value of K 2 . Away from the light 
cone, thermal corrections of order 1/mT appear in I{K). 



B Relations between some integrals 

In this appendix, we establish some useful relations between the five integrals 
defined in Eq. ©. Let us start by the study of J AAB '■ this integral satisfies the 
2x2 linear system 

' i W — T 

" AAB AB (39) 

B 1^ — J 

P AAB AA 1 

the resolution of which gives the two components of J AAB as functions of J AA 



and J AB . 



ftP J — AP J 



In order to obtain J AABB it is convenient to define first the second rank tensor 
aabb ~ J 2w 2loi (L+ ■ A) 2 (L+ ■ B) 2 [ ' 

— oo 

from which we can obtain J AABB as 9^vJ A u ABB - The three independent compo- 
nents of this symmetric tensor can be obtained via the resolution of the following 
3x3 linear system 

'^B v J^ abb =J AB 
A p A ^Zbb = Jbb (42) 



Bjj,B v J AABB ~ J AA 1 



which finally gives 



2(A-B)J A B -A 2 J A A -B 2 J j 
(AAB) 

The same method can be applied to J AABC , which gives 



T V J AB A A BB 1 AO\ 

AABB — (An D\2 ' V ' 



j _ A 2 J AA + B 2 J AB + C 2 J AC 

aabo (AAB) 2 ' [ ' 
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under the assumption that A + B + C = 0. We also need the following relation 

Jab + Jbc +J OA =Q> ( 4 5) 

which is valid when A + B + C = 0. 
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